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We investigate how far a new physics scenario affecting primarily the third generation fermions can 
ameliorate the tension between B-decay observables and Standard Model expectations. Adopting 
a model-independent approach, we find that among the three observables that show signs of such 
a tension, viz. the branching fractions for tv, Bd — > D{D*)tv, and the like-sign dimuon 

anomaly in neutral B decays, the first two can be explained adequately, while there is only a marginal 
' improvement for the third. As a spin-off, it is shown that one can also accommodate a change in the 

, branching fraction of the Higgs boson to a r lepton pair from the SM expectation, if such a change 

■ is estabhshed in future data. 
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I. INTRODUCTION 



While the purported discovery of the Higgs boson at the Large Hadron Collider (LHC) (IJ, |2j] seems to vindicate the 
(— I ] Standard Model (SM), there are enough reasons to believe that the latter is but an effective theory, valid only up to 
I— !■ a certain energy scale, with a more complete theory lurking beyond that. One of the major reasons for such a belief 
is the fine-tuning problem associated with such an elementary scalar; also of considerable import are issues such as 
^ ■ the existence of the dark matter, or the baryon asymmetry in the universe. 

' This acts as a strong motivation to look for signals, both direct and indirect, of such a new theory. While direct 
, signals very often involve production of new particles, the indirect signals will, most probably, be manifested as 
■ modifications to SM observables by new effective operators, or even the same operators as in the SM, but with 
] modified Wilson coefficients. 

B meson observables have constituted a favourite hunting ground for indirect signals. Over the years, several 
experiments, including B-factories, Tevatron, and even the LHC, have reported observables that are not in good 
. agreement with the SM. While the tension is not so overwhelming as to claim unquestioned evidence of New Physics 
(NP), the pattern is interesting. Here, one must remember to tackle the theoretical uncertainties carefully; some of 
the discrepancies, like the longitudinal polarization anomaly in the decay of a B meson to two vector mesons, vanished 
because of a more careful reappraisal of the SM effects. 

Let us begin by considering a few observables which are not in full conformity with the SM expectations; 



• the large branching ratio oi B ^ D{D*)tv, with a combined tension of 3.4(T Q; 

• the large branching ratio of r+i^, with a tension of 1.6cr [4| ^; 

• the like-sign dimuon asymmetry, with a tension of 3.9tT 

It is interesting to note that the first two involve a r lepton in the final state. This motivates us to ask if there exists 
one or more new effective operators involving the b quark and the r lepton. Such a possibility was raised in Ref. , 
and further investigated in Refs. One might feel tempted to add to this list a hint of another anomaly: the 

branching ratio oi H ^ t^t^ seems to be a bit on the lower side than that expected in the SM 
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^ The tension has come down very recently; it was about 2.8<t before the publication of the latest Belle result [^. 
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The lowest dimensional operators of interest can, generieally, be expressed as (bTAs) (tTbt) where Ta,b are 
appropriate combinations of the Dirac matrices. As shown in Ref. 0], this four-fermion operator is relatively uncon- 
strained. This leads to a new contribution to the Bg-Bg mixing amplitude, with a nonzero absorptive part: cutting 
the intermediate r propagators can yield an on-shell r-pair. Thus, one has a new contribution to AFs, the difference 
in the widths of the two B^ mass eigenstates, which, in turn, ameliorates the apparent discrepancy in the like-sign 
dimuon asymmetry. However, the strength of any such operator is ultimately constrained by the mass difference A.Mg 
of the Bg mass eigenstates. 

Considering the fact that there is hardly any tension in the data involving electrons or muons in the final state 
(except the dimuon anomaly, to which we will come shortly), one might feel tempted to invoke one or more effective 
operators involving only third generation quarks and leptons. While such an effective operator based study was 
undertaken in Refs. the constraints on AMg were not correlated with those coming from the decay width 

difference ATs] they were assumed to be independent numbers. The authors of Ref. discussed the effectiveness 
of AMg as a possible constraint on the parameter space. 

In this paper, we adopt a different approach. To begin with, we posit a single effective operator involving a 
third-generation quark current and a third-generation lepton current. As it involves only third-generation fields, 
the constraints on the Lorentz structure for the same, or on the magnitude of the corresponding Wilson coefficient is 
relatively weak. For example, just restricting the new couplings to the perturbativc regime ensures that T(IS') — > t^t^ 
does not receive a significant contribution over and above the SM amplitude, which is electromagnetic in nature. We 
will, however, not venture to discuss any particular models that might predict such an interaction, and adopt, instead, 
a bottom-up approach. 

A theory of flavour would, generieally, dictate that such an operator would be written in the weak basis. On the 
breaking of the electroweak symmetry, the quark fields would need to be rotated to the mass basis. This leads to a 
plethora of new operators, related to the original through the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. 
Included, amongst others, are those leading to 6 — ?> sT+r~ and b svv (the lepton fields are not rotated, so we 
will always write v for Vr). We will see that, even with a moderate NP scale ^ 1 TeV, the constraints from AM^ 
are so strong for vector and axial-vector ({V,A}) operators that the effects on AF^ and Bg — )■ r+r" are bound to 
be unobservably small. Thus, the explanation of the dimuon anomaly must lie elsewhere, while this scheme can 
successfully explain the charged-current decays. The outlook is better if the effective operators are of scalar and 
pseudoscalar ({5, P}) variety. Indeed, such a scenario predicts a rather strong enhancement of the branching fraction 
of Be — >■ TV over its SM prediction. As for the tensor current operators, the corresponding Wilson coefficients are 
severely constrained from radiative decays like b — > 57. and so we will not consider them any further. 

The rest of the paper is arranged as follows. In the next section, we will briefly go through the existing data. In 
Section III, we will discuss the new operators; flrst, the {5*, P} type, and then the {V, A} type. In Section IV, we show 
how these operators may help in bringing down the tension with the SM. We summarize and conclude in Section V. 



II. EXISTING CONSTRAINTS 



A. B ^ D{D*)tu 

The importance of studying the B — > D{D*)tv modes for a possible signal of new physics has already been pointed 
out in the literature [l^. The BaBar Collaboration 0] measured the branching fractions for these two modes, and 
they are above the SM predictions. They are also not consistent with a type II two-Higgs doublet model (such as the 
minimal supersymmetric extension of the SM) . The implications of the data as a possible hint of physics beyond the 
SM have been studied in [T2j . 

It is particularly useful to consider the ratios R{D) and R{D*), defined as 

^ ' Br(P ^ P>(*)£i/) ^ ' 

as these are largely free of the uncertainties — e.g., those in the form factors — that exclusive modes are often prey to. 
The SM predictions are 

P(i:>) = 0.297 ±0.017, P(i:»*) = 0.252 ±0.003, (2) 
while the BaBar Collaboration quotes 0] 



^(1?) = 0.440 ± 0.058 ±0.042, P(L»*) 0.332 ± 0.024 ± 0.018 . (3) 
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It should be noted that a recent calculation [l^ in unqucnchcd lattice QCD gives, in the SM, R{D) = 0.316 ± 0.012 ± 
0.007. This is consistent with the earlier SM prediction, but cannot explain the tension with the data^. 
Using Eqs. ([2]) and ([3]), and adding all errors in quadrature, we get 

= 1.481 X (1 ± 0.173) , = 1.317 x (1 ± 0.091) . (4) 

B. B ^ TV and — ^ tv 
The partial decay width B — > tv, in the SM, is given by 



,2 
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The world average is 01 

Br(B ^ Ti/) = (11.5 ±2.3) X 10"^ (6) 



while the theoretical prediction is 



Br(B ^ Tiy)su = (7.57^^61) x W , (7) 



which gives a tension at the level of 1.6a 1J|. The theoretical uncertainty comes from those in the decay constant 
/b and the CKM matrix element Vub- While the discrepancy has eased considerably, from 2.8a to 1.6a, after the 
publication of the new Belle result there is still a non- negligible tension between the value of \Vub\ determined from 
this decay, and that determined indirectly from the sides of the unitarity triangle, or an average of direct inclusive 
{B — )■ Xuiiy) and exclusive {B Triv) measurements. 

The discrepancy has led to several attempts in the literature to explain this as a possible NP signal. However, the 
explanations based on the existence of only a charged Higgs boson of type-II are ruled out at 95% confidence limit 
from the combined analysis of processes like b — >■ S7, Z — >■ bb, B^ — B^ mixing [l3l; the goodness-of-fit is hardly better 
than the fit with the SM alone. Models like R-parity violating supersymmetry fare much better and give a satisfactory 
explanation of the excess [l^. Note that all those analyses were performed using the pre-2012 data. 

A similar expression as in Eq. ([5]) holds for B^ — >■ tv. For numerical evaluation, one might use 

/s, = (395 ± 15) MeV , t^^ = 0.458 ± 0.030 ps . (8) 

C. Bs—Bs mixing: AAIs, ATs, Ps and 0s 

While there are no apparent tensions in this sector at present, the data, as we will soon see, acts as a very tight 
constraint on NP operators. The mass splitting between two Bg mass eigenstates, AMs ~ 2|A/i2s|, is extremely 
well- measured [l^, namely 

AMs = 17.719 ± 0.043 ps^^ . (9) 

This agrees very well with the SM expectation , viz. 

AMs (SM) = (17.3 ± 2.6) ps"^ , (10) 

and acts as a very tight constraint on NP models. There are two relevant phases in the Bg-Bs system. The first one, 
the mixing phase, is defined as 



^ The Belle collaboration measurements fisll viz. 

R(D'\ - 70+019+011 R(D*^ - 47+0-11+0.06 

while being even further away from the SM expectations, nonetheless are consistent with these as well as the BaBar results. This 
agreement, though, is but a consequence of the large error margins. 
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while the second one, responsible for semileptonic asymmetries, is given by 

0. = arg(-^). (12) 

The SM predictions [Hi are 

(j), = 0.0041 ± 0.0007 , ~2Ps = -0.038 ± 0.002 . (13) 

The experimental numbers are 

-2/3, = -0.040™ (14) 

from direct determination, and 

- 2/?, = -omest'oZit (15) 

from an indirect global fit [l^. We will use the former number. 

Note that Bg lifetime is rather ill-defined, as the two mass eigcnstates have a significant lifetime difference, namely 

Ti3^^ = 1.408 ±0.017 ps , Ti3,„ = 1.626 ± 0.023 ps. (16) 

Averaging over the two, we have 

2 

TB, (average) = — = 1.509 ± 0.012 ps , (17) 

f L + f if 

which should be compared with tb^ ~ 1.519 ± 0.007 ps. Thus, 

— = 0.993 ±0.009 , (18) 

while the SM expectation for this ratio lies between 0.99 and 1.01 [T9| . 
The width difference AF, is given by 

Ar, =2|ri2,|cos(0,). (19) 

While the SM predicts AFj > 0, there was a sign ambiguity earlier in its determination. Recently LHCb, from the 
decay Bg — )■ J /i^K'^ K~ , found AF, > with a 4.7a confidence level [2l|. The experimental number, an average over 
various measurements |19l |. 

Ar, = 0.095 ±0.014 ps-i (20) 

is to be compared with the SM prediction 

Ar, (SM) = 0.087 ± 0.021 ps^^ . (21) 

D. The like-sign dimuon asymmetry 

The like-sign dimuon asymmetry, defined as 

, _ N{fi+^i+)-N{ti-fi-) 
^«'-iV(M+A^+) + 7V(M-M-)' ^ ' 

and measured with 9.0 fb~^ of data at the D0 Collaboration is @ 

All = (-7.87 ± 1.96) X IQ-^ . (23) 
This can be expressed as individual flavour-specific (fs) semileptonic asymmetries coming from Bd and Bgi 

All = (0-595 ± 0.022) aj^ + (0.405 T 0.022) a^^ , (24) 
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where the numbers in the parentheses are the production fractions for Bd and Bg. The SM expectations are 

a^, = (-4.1±0.6) X 10"^ a}, = (1.9±0.3) X 10"^ (25) 

which gives the SM prediction 

(A^;)sA/ = (-2.4±0.4) X 10-^ (26) 

Comparing Eqs. (j23p and (j26p . one finds a 3.9cr discrepancy between theoretical prediction and experiment, a^^ has 
aheady been measured by BaBar and Behe; the average [l^ 

= (-3.3 ±3.3) X 10-3 (27) 

is consistent with the SM. This gives an indirect prediction for a^^, viz. 

a^^ (-1.52 ± 1.04) X 10-^ (28) 

where the error has been symmetrized. We have neglected the correlation between a^^ and a^^, but have taken 
the uncertainties in the production fractions into account. Recently, the D0 Collaboration directly measured a^j, = 
(-1.08±0.72 (stat)±0.17 (syst)) x 10"^ [12] which is also consistent with the SM expectation. The HFAG collaboration 
averages over several direct measurements of a^^ and quotes [l9j 

a}, = (-1.05 ± 0.64) X lO^^ (29) 

but this has a nonzero correlation with a^^ . 
This gives a weak constraint on 0^: 



tan0, =a},^. (30) 



AM, 
AT 

If there is some NP contributing to both Mi2s and Ti2s, one can parametrize the NP contribution as 

A/i2 = i\/f/^ + i\/i^^ = Aff2*^i?M exp(^0M), 
i 12 — i 12 + i 12 = i 12 exp(_i0rj , 

resulting in Q 

4>s = + 4>M - <t>T . (32) 

Thus, there are two ways to have a large a^^; either a large contribution to APs or a large - 7r/2. But 0f^^ w 
and (pi^i = —2j3s is known to be small, so a large (j)s necessarily warrants a large (pr, and hence a large contribution 

to Pi2. 

Taking all the existing constraints into account, it was shown [3,B[ili that b — > st^t^ is a viable option to generate 
a large P12. However, such a new channel decreases the lifetime of Bg compared to Bd] moreover, one does not expect 
Br( i?s — >■ T+T") to be more than 3-3.5% [ll|. The inclusive mode B{Bd XsTt) is constrained to be less than 5% 
[ig, while BaBar gives a 90% limit ^ 

B{B+ ^ i^+r+r-)|,2>i4.23 GcV^ < 3.3 x lO'^ . (33) 

However, not all Lorentz structures that contribute to a new absorptive part in Bg-Bs mixing contribute simultane- 
ously to Bs T^T^ or B^ K^T^T^ . At the same time, there can be significant long-distance effects in Bg-Bs 
mixing, through meson loops, and they can have a non- negligible contribution in APs 25,[26|. While we will discuss 
these issues in detail later, the crucial point to note is that the NP contributing to P12 should, in general, contribute 
to M12 also, and the mass difference AAIg is so tightly constrained that this leaves only a very small room for any 
NP. 



E. Di-tau suppression 

Looking for the SM Higgs in the H ^ tt mode, the CMS collaboration has failed to see [l[ an unambiguous excess 
over the background. Indeed, for the preferred mass of mn = 125 GeV (corresponding to the much-touted diphoton 
and four-lepton excesses), it is only able to impose a 95% C.L. upper limit on the ditau excess. With the result being 
similar for the ATLAS collaboration Q as well, the pp ^ H ^ tt cross-section is, in fact, consistent with zero, viz.. 
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Observable 


SM 


Expt 


-2/3, 

tbJtb^ 
AMs « 2|Mi2,| 

Ar, f» 2|ri2s|cos0, 

a% = all 


-0.038 ±0.002 

0.99 - 1.01 
(17.3 ±2.6) ps-^ 
(0.087 ±0.021) ps^^ 
(-2.4 ±0.4) X lO""* 
(-4.1 ±0.6) X lO""* 
(1.9 ±0.3) X 10"^ 


_n 040+°'''"' 

U.U4U_Q Qgg 

1.509 ±0.012 ps 
0.993 ± 0.009 
(17.719 ±0.043) ps-^ 
(0.095 ± 0.014) ps^^ 
(-7.87 ± 1.96) X lO"-' 
(-3.3 ±3.3) X 10"^ 
(-1.05 ±0.64) X 10"^ 


R{D) 
R{D') 
Br(B+ ^ r+u) 


0.297 ±0.017 
0.252 ± 0.003 
(7.57«:«,?) X 10-^ 


0.440 ± 0.072 
0.332 ± 0.030 
(11.5 ± 2.3) X 10"^ 



TABLE I: Inputs used in our analysis. For details, see text. 



F. Numbers used for the analysis 

Apart from the numbers shown in the previous subsections, a summary of which is given in Table 1, we also use 
the following for our analysis: 

mB+ = ms, = 5.279 GcV , tb+ = 1.641 ps , rs, = 1.519 ps , tob, = 5.367 GeV , (34) 

and 

\Vtd\ = (8.67t° j?) X 10-3 , \Vts\ - 0.0404t° °°ii , |K,| = OMUt'oZl , 

= (3.49 ± 0.13) X 10-3, 7 « arg(l/*6) - 77° . (35) 

Note that \Vts\ is measured from Bg-Bg mixing, but if we talk about new physics in the mixing and hence AAfg, 
we should, instead, use \Vts\ as determined from the unitarity constraints. The central value as determined from the 
unitarity is 0.0404; purely from AMg measurement, this comes out to be 0.0429 ± 0.0026. The error margin in 7 
is not important for our analysis, so we take the central value [l^. For the evaluation of AMg, we have used the 
unquenched lattice value 



/b.a/BbT^ 248± 15 MeV. (36) 



III. THE NEW EFFECTIVE OPERATORS 



Let us, now, consider a set of possible operators involving third-generation fermions, satisfying both Lorentz and 
gauge invariance. These might be of the form 

Os = AiQs^dsRXesRLsL) + B{Q^lUzr){^:,rL'Ii^) ± h.c. (37) 

with L'3 = 4(72^3, or 

Ov^C [03i7^T,g3L] [X3L7M^ai3L] , (38) 

where Q3, L3, U3, o?3, and 63 stand for the doublet quark, doublet lepton, singlet up- type, singlet down-type, and 
singlet charged lepton of the third generation respectively. In view of the experimental measurements that we seek 
to address, we limit ourselves, here, to only those operators that admit charged-current interactions. Furthermore, 
we do not consider tensor operators as their Wilson coefficients are very tightly constrained from radiative decays. In 
terms of component fields, we can write the scalar-pseudoscalar operators as 

Os = A[(6Pfl6)(rPLr)±(<Pfl6)(TPLt^)±(foPL6)(TPflr)±(5Pit)(77P;^r)] 

+B [(tPRt){TPLT) - {bPLt){VPRT) ± {tPLt){TPRT) - (tPRb){TPLu)\ , 
= A [ibPRb)(TPLT) + h.c] ± B [{tPRt){TPLT) + h.c] 

+ {A - B) [(t.PRb)(TPLiy) + {bPLt)(vPRT)] 
= I- A [(66) (tt) — (6756) (r75r)] + similar terms . (39) 
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Eq. (|39l) shows that the neutral current operators have a coefficient different from that for the charged current 
operators. In fact, there are two neutral current operators now, one involving scalar currents and the other involving 
pseudoscalar currents. While we will discuss later the consequences of such an operator structure, note that the r 
Yukawa coupling can, in principle, be significantly modified by a top loop. Corrections to the Yukawa couplings of 
other third generation fermions are negligible. 
In a similar vein, 

Ov - C[{h^'PLt){V-1^PLT) + {trPLh){T-1^PLv) 

+ \ (pl^Phb - h''PLt){Tj^PLT - Vj^Pli^)] . (40) 

To make explicit the higher-dimensional nature of the couplings, we denote 

A^a/A\ B^b/A\ C^c/A\ (41) 

where a, b, and c are dimensionless couplings. 

If Eqs. ([5^ and/or PO)) are all we have, the phenomenology is straightforward, and only a subset of that we would 
consider below. One might think that T — > tt will put a tight constraint on the coefficients, but, in actuality, that 
constraint is far too weak. The reason is that the SM decay is an electromagnetic one, and the width is given by [l^ 

rT(i5)^« = 4 a2 q2 ^,^-2 |^(o)|2 ^ 2x) yr^, (42) 
where x = /M^ and i?(0) is the radial part of the non-relativistic wave function at the origin. 



A. The rotation 



Assuming that the operators in question have arisen on account of some flavour physics operative at scales higher 
than the weak scale, we now put forward the Ansatz that the fields in Eqs. (p9| and (|40|) are in the weak basis, 
and should be rotated to the stationary or mass basis. Let us, for simplicity, assume that right-chiral fields are not 
rotated, and for the left-chiral fields, 

bnik xid + X2S + X3b , t^k ~> yiu + y2C + Hat , (43) 

where the right-hand side fields are in the mass basis. 

If U and V matrices are responsible for the rotation of T3 = -|-l/2 and T3 = —1/2 fields from the weak basis to the 
mass basis, so that the CKM matrix V = U^V, one notes that {xi,X2,X3) and (j/i, 2/2,2/3) are just the third rows of 
2? and U respectively. If we assume the rotation matrices to be almost diagonal, the only constraint is 

y*3X3 « Vtb . (44) 

As for other combinations, we can, at most, use order-of-magnitude arguments to yield 

y^xi ~ Vtd , 2/3X2 ~ Vts , 2/22:3 ^ Vcb , ylx3 Vub , (45) 

although there can be significant deviations. 

IV. THE OBSERVABLES 
A. Leptonic and semileptonic decay channels: t^J^j Be f^i^, B — > D{D*)tu 

The relevant new operators are 

Os D j^{a~h)[{{ylu + y*2c}PRb) (rP^^^) + h.c] , 

Ov D ■^cx3[i{ylu + y;c})PLb) (rPiJ/) + h.c] , (46) 
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FIG. 1: The dependence of Br(_B tv) on the magnitude of the new physics couphngs |cj/*2;3| (thick black bands) and 
I (a — h)yi\ (thin purple/light grey bands). In each case, the different bands from left to right correspond to differing values 
of the phase of the new physics coupling, namely, — 7r/2,0, 7r/4 and 7r/2 respectively and the thickness of the individual band 
reflects the errors in |Kib| and fs at the Icr level. The experimental data on Br(i3 — >■ rv) at Icr (red/dark grey lines) and 2cr 
(blue/light grey lines) intervals are shown as horizontal bands. 



and their effect on tlie amplitudes of interest can be obtained by simple replacements in the corresponding SM 
expressions, namely, 



7f Vf 4A2^i^3 forOv, 



(47) 



where the latter follows from 



{Q\u{l-j5)b\B-) = -ifB (48) 

nib + niu 

For the decay, one has to make the following substitutions: {u, Vub, 2/i , , m„ , Jb} {c, Vcb, 2/2, ^ris^ , rtic, /bA- 

Taking only one set of new physics couplings, ca^ay* or (a — b)yl, to be non-zero, in Fig. [T]we show the variation 
of Bt{B — >■ Ti/) with the magnitude of the coupling. We have set the scale of the new physics A = 1 TeV, and used 
4>3 = J [13 • Understandably, the phase of the NP coupling plays a significant role with positive values allowing for 
destructive interference with the SM amplitude. This results in the different bands (one for each representative value 
of the phase). The width of the bands is due to the uncertainty in Vub and to a lesser extent, that in /b- The two 
horizontal bands correspond to la (red/dark grey lines) and 2a (blue/light grey lines) intervals of experimental data 
on Br(i3 — > Tf). Their intersection with the NP bands determine the allowed ranges for the couplings. Note that if 
\yi\ is indeed OdT^bl), then |c?/*a;3| > 0(0. 1) would indicate a significant departure from the expectations in naive 
dimensional analysis. 

For B D{D*)tv and Be — > tv, the SM effective Lagrangian is 

C^fj = in'PLb) {t^.Pli^) . (49) 

The two NP operators Os and Oy modify this to 

Cmp = [(1 + g') icrPLb) {T-f,PLv) + gn, (cPRb) (tPlu)] , (50) 

where the g' (gR) term emanates from Oy (Cs), namely 

^ g = Cy2X3 , ^ gB. = [A- B)y2 . (51) 

If (?' 7^ but 5_R = 0, one can write 

R{D) = RsM{D)\l+g'\^ (52) 
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FIG. 2: Allowed values of the magnitude and phase of the new physics coupling cyiXs. The black, green (light grey), and red 
(dark grey) regions denote the la, 2a, and 3a consistency bands respectively, where both theoretical and experimental errors 
are taken into account and added in quadrature. 



and a similar equation for R{D*), assuming that the new interaction does not contribute to the electron or muon 
channel. On the other hand, if gj^ ^ but g' = 0, we get 

R(D) = RsM{D)[l + 1.5Rc{gR) + \gRf] , 
R{D*) ^ i?sM(i?*) [1 + 0.12 Re(gfl)+ 0.05 . (53) 

The same couplings also contribute to the leptonic decay Be — > tv, and depending on the phase of the coupling, can 
increase or decrease the branching ratio. 

In Figure[2l we show the allowed values of the coupling cy2Xz, with A = 1 TeV, at different confidence levels. The 
intervals are calculated with individual error margins and not with a combined X^-fit. Corresponding to these three 
cases, the branching ratio of Be — > tv are 



Br(Bc ^ Ti^) e [2.05 - 2.40]% (Ict) , [1.80 - 2.60]% (2cr) , [1.60- 2.80]% (3ct) , 



(54) 



which should be compared with the SM value of 1.68%. 

For the {S,P} couplings, there is no region in the parameter space compatible with both R{D) and R{D*) at la 
or 2cr level. This is because of the small contribution of gn to R{D*). Only at 3cr, does one get an allowed region in 
the parameter space. However, Br(i?c tv) can be quite large because of the chiral enhancement. If we assume, as 
a conservative estimate, Bi{Bc — > tv) < 10%, this translates to 



I (a-6)y* I < 1.05 



(55) 



for A 



1 TeV. Note that the limits one obtains from neutral current mediated processes, like Bg-Bg mixing or 
'T~ [ll[, are not valid in these cases as the couplings are different. 



B. Bs T+T- 



The term from Os that we would be interested in is 

A [{{Xld + Xls + xlb}Plib){TPLT) + {bPL{xid + X2S + Xzh}){TPRT)\ 

and the relevant operator is 

Axl{-sPRb){TPLT) = ]^AX*2 [isPRb){TT) - isPRb)iTj5T)] . 



(56) 



(57) 
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This gives 



where 



A = 
B = 

c = 



2567r3 



Anii 



mi 



Am 



1 ^ \A + B 



(58) 



C 



mfc + ms 

c 



2m-r 



mb + ms mg 



\pl An 



(59) 



where Cio is the Wilson coefScient for the corresponding SM operator, and is too small to be of any consequence. If 
we take the upper limit of Br(_Bs — > r+r^) to be 3.5%, we get a bound, depending on definition used for mb. For 
example, 



.;i< 1.52 (1.34) ( 



3.5% 



1/2 



(60) 



for mb = 771^°'*^ = 4.8 GeV {mt = mb{mb) = A.2 GeV). This agrees with Ref. [ij. Thus, potentially, a can be large. 
For the {V, A} couplings coming from Ov, one has 



Bl{Bs^T^T ) = — Jl — 

327r \/ ml 



1 c , 



Note that this is further suppressed by a factor of m^/m^ , which results in a weaker bound: 



3.5% 



(61) 



(62) 



Width difference ATs 



Following Ref. [ll|, let us quote the relevant expressions for the width difference AFg 

F? 



where 



Os 
Ov 



12, NP 
^12, NP ~ ^ 



N 



{Q§) 



3Nxy/l~Ax{Q§)Cl, 

{l + (l-x)^/^^} (Q^) 

+ {l + (l + 2x)VT^} (Qf) 

2 / 2 
-7 \VtsVtb) , 

67rmSa 
5 

" 12 



Cs 
Cv 



V y 



1 



y2_ 

AGf V*,Vtb 2A2 
V2 1 c 
AG^Wv^b^^ 



(63) 



(64) 
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Note that the second equation of (|63|) has to be augmented by the inclusion of the Vt loop, which can be obtained from 
the corresponding r contribution by putting x = 0. For numerical evaluation, we use the lattice values = 0.231 
GcV, Bg = 1.3, By = 0.84. From B — > A'rr, there is a (scale-independent) bound, namely Cy < 0.8, which translates 
to 

^cx^xa < 1.05 (^p^) . (65) 



D. AMs and the mixing phase (/)m 



The aim of this subsection is to show how and why the constraints coming from AM., measurements are so restrictive 
in nature. Here, we will start from Oy- There can be two possible diagrams, one with the r lepton and the other 
with the neutrino. As the amplitude is not chirality suppressed, both the diagrams contribute equally. The exact 
amplitude cannot be calculated unless the know about the ultraviolet completion of the effective theory. If we use 
a cut-off regularization, the leading term, which is divergent, should match with the leading term of the full theory. 
The leading term of the loop amplitude is quadratically divergent, so we can safely neglect the sub leading terms. 

The relevant part of the effective operator is 

Ov D lc{h''PLb){T-ff,PLT) ^ ^Cx*X3{s^''PLb)iT-f^PLT) . (66) 



This gives rise to a mixing amplitude 



iMi2s = ( -^cx^xs ) (Ol)(^^2) x 16 



^icx*x,fA-'{Oi) 
^ {cx^x^f K-'^vbMbJIBb,, . (67) 



2 



where we have used 



2 



r2 = ^, (Oi) = -t/bA/b/^Bb , (68) 

r2 being the leading term of the loop amplitude, and Oi = [sa7'^(l — 75)6^] [s^7^(l — 75)6/3], a and /3 being colour 
indices. The factor of 16 can be understood in the following way: there is another crossed box, which, in an effective 
theory, is something like a t-channcl amplitude. This gives a factor of 2. The initial meson can pick up a 6 from Oi 
in two ways, and an s in two ways, so the symmetry factor is 4. The neutrino mediated -best amplitude gives another 
factor of 2. 

Comparing with iMf^J , we find 

A{cx*2X3)^A- ' 



(69) 



where xt = m^/rn^^ and So{xt) is the relevant Inami-Lim function. The SM amplitude is GIM suppressed whereas 
there is no such suppression for the NP amplitude, and thus Eq. (j69p puts a fairly tight constraint on 0X2X3. If we 
want the latter to be large, the phase must be opposite to that of the SM amplitude, so that there is a destructive 
interference: Af^^ ^ ~2M^^^ . Taking the errors on the AMg prediction in the SM and the measurement of the 
same quantity both at 2a, we get 

\cx*2X3\ < 0.048 f ) . (70) 



1 TcV, 

The allowed region is shown in Fig. [31 For A = 1 TcV, the limits on a^^ are 

- 6.3(-11.2) X lO"'^ < a}^ < 2.2(6.9) x 10"'' (71) 

at 1(2)ct. Thus, by themselves, such operators are unable to explain the dimuon anomaly, and the explanation must 
lie somewhere else. The maximum value of Br(i?s — r+r") is about 3 x 10""'%. 
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Onp/JI 
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FIG. 3: The allowed magnitude and phase of the coupling 0X2X3 from the measurement of AMs and 2/3s, the effective mixing 
phase from the box amplitude. The light (green) shaded area between the dashed (blue) curves is allowed by AMs measurement, 
whereas the area between the solid (purple) curves is allowed by the data on 2/3s. Thus, only the two patches of dark (pink) 
shaded area is finally allowed. For the left (right) plot, the experimental errors are taken at 1(2)(t level. 





Hiq) 

t(p) t(p) 

t{k + q) 

FIG. 4: Typical one-loop corrections to (left) the Htt vertex and (right) r self energy. 



The situation is marginally better for a chiral coupling in the scalar sector; i.e., either S — P 01 S + P. For such 
operators, the leading term in the Bg-Bg mixing amplitude is proportional to m^logA^, and there is no effective 
constraint from AA/s and mixing phase. However, the major constraint comes from Br(i?s — > r+r^), and also 
partially from AF^. Wc find 



\a)g\ < 6 X 10" 



(72) 



for Br(i3s — > r+r ) < 4% and taking all the errors at 2a level. Thus, none of these schemes are enough to explain 
the dimuon anomaly completely. 



E. H 



Let us begin by parametrizing the tree-level Higgs-tau coupling by 



{Htt) 



h-r T T H 



(73) 



As in any quantum theory, this interaction Lagrangian receives quantum corrections. We neglect here all the SM 
corrections and concentrate solely on that wrought by the four-Fermi operators. The scalar-pseudoscalar operators 
give rise to an effective interaction of the form 



^Re(y3) [{tt)(rT) - (^TsOC^Tst)] + ^Mys) [(tt)iri5r) ~ Mirr)] 



2A2 



(74) 
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Each of these terms generates, at one-loop, a two-point diagram contributing to the effective Htt couphng (see 
Fig. |4]). The said diagrams are manifestly quadratically divergent and need to be regularized. Given that our basic 
theory is only an effective one, we may use a momentum cut-off regularization scheme, to yield the following correction 
to the Lagrangian of Eq. ([75)) : 

-^^i-ioop = ^ [Re(2/3) rr + lni{y:,) t^^t] H + ---, (75) 

where the ellipses denote sub leading terms. It is natural to consider Acutoff ~ A, for the two are expected to be 
similar. The appearance of a divergent correction to the pseudoscalar coupling (one that did not exist at the tree 
level) might seem disconcerting at first. However, it should be realised that we are dealing with a nonrenormalizable 
theory and the existence of such a divergence only reflects the fact that a large correction to Htt is not prevented by 
the symmetries of the theory extant on admitting the general four-fermion interaction. On inclusion of the ultraviolet 
completion, such divergences would disappear identically. Clearly the two Lorentz structures contribute incoherently 
to T{H — >■ rr). Formally though, the contribution of the scalar coupling correction may be larger as it can interfere 
with the SM amplitude, and thus can enter at an earlier order in the perturbation theory. For simplicity, though, let 
us assume that is real. Then, we can parametrize the effective Htt vertex, upto one-loop by 

4^'^ - hr{l+OrrH+--- 

Similarly, several other Yukawa couplings also receive corrections, but these are suppressed on account of the particular 
structure of NP. For example, the bottom quark Yukawa coupling receives a correction on account of a tau-loop, and 
this change can be expressed as hb ^ hb + bhrVs/STT^. 

It might be argued that our calculation of ^ is somewhat naive, and it is indeed true. However, an exact calculation 
necessitates a knowledge of the ultraviolet completion of the theory, and, an a sense, goes against the spirit of an 
effective theory. Nonetheless, ^ does encapsulate the leading correction, and in Fig. [Sj we show the variation of the 
branching fraction oi H ^ t^t~ as a function of the real variable ^ for mn = 125 GeV. From Fig. [5l it is very clear 
that even a moderate value of ^ ~ —0.3 is enough to give a 50% suppression in the BR(iJ — > t+t^) . On the other 
hand, from the observed upper limit of gr^tt I o'SM ~ 1-1 from the LHC [l|, we get an upper limit of ^ « 0.05. 
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For the {V, A} current, once again, both scalar and pscudoscalar couphngs appear at one-loop. However, the loop 
is convergent as the current structure demands that extra powers of the fermion masses must be picked up. Hence, 
the corresponding corrections arc too small to be of any consequence. 

F. Anomalous top decays 

One might ask whether the new couplings will lead to observable rates for FCNC top decays, e.g. t — > ct+t^. 
Unfortunately though, even if wc use values of the couplings 61/2 (for scalar operators) or (for vector operators) 
significantly larger than what wc need to explain the anomalies under investigation, the rates for this decay are still 
much smaller than the LHC sensitivity limits. For example, we might naively use the limit on the branching ratio of 
t ^ cZ namely 

Br(< ^ cZ) < 0.24%, (77) 

along with the measurement [30| of the decay width of the top 

Tt = 2.00tg;g GeV (78) 

to yield the very weak limit of A^0.5to4 when the couplings are only restricted to be perturbative. Furthermore, 
even the use of Ea. (|77p is over-optimistic, for the CMS limits have been derived requiring that the Z-mass can be 
reconstructed from its decay products. In the current case, this does not apply and the signal to noise ratio is lower 
than that assumed to obtain Eq. ([77| . In other words, the actual limit is much weaker than that quoted above. 

V. CONCLUSIONS 

In this paper, we have investigated the possible implications of a scenario that involves some new interactions 
involving the third generation fields. Any model that treats the third generation differently from the first two 
generations may lead to such a scenario. Without attempting to prescribe an ultraviolet-complete theory, we rather 
consider an effective theory valid below some cutoff scale A, above which the full theory takes over. A possible 
motivation for such a scenario is the fact that there are excesses over the SM predictions for the charged current 
B-decays, namely B D{D*)tv and i?+ tv, while the predictions for the processes involving the first two 
generations of leptons do not show any tension with the data. 

In the effective theory, there can be several four-fermion operators involving the third generation fields and several 
possible choices for the Lorentz structures of the currents. With the Wilson coefficients for tensor operators being 
severely constrained by the data on radiative decays, we preclude these from our discussions. With A being larger 
than the clcctroweak scale, it is quite likely that such four-fermion operators in the effective Lagrangian should be 
written in the weak basis, and for reasons of economy, we consider only one such operator at a time. Rotating the 
fields to the mass basis generates new operators involving first and second generation quark fields, albeit suppressed 
by the corresponding entries of the quark mixing matrix. 

Once we have a set of such operators, wc study their implications on several B-dccay obscrvablcs. In particular, 
wc show that the apparent excesses in the B-dccay channels mentioned above can be accommodated satisfactorily in 
this scenario; complementary observables lead to nontrivial constraints on the model parameters. The vector-axial 
vector operators successfully explain the excesses in both B — > Dtv and B — )• D*tv channels, apart from leading to a 
sizable enhancement to the Be — > tv branching ratio as a testable prediction. The scalar-pseudoscalar couplings are 
not that successful in explaining both the excesses, but there is a definite improvement over the SM predictions. The 
excess in the channel B'^ — >■ tv can have a satisfactory explanation too, although the tension is no longer worrying. 

The operators leading to Bd-Bd and B^-Bs mixing are more constrained. They have identical Lorentz structures 
as those discussed before, but with different quark fields and different Wilson coefficients. While these coefficients 
are constrained from the measured mass differences AM^ and AMg, the restrictions are not strong enough to rule 
out any observable enhancement in the Bg — > t^t~ channel, which should be investigated more carefully as one of 
the best windows to new physics. Unfortunately though, the anomalously large dimuon asymmetry receives only a 
marginal improvement over the SM prediction. One might need other operators to explain this, but it is not easy 
given the tight constraints from AAfg measurements. 

One thing that still remains unobservably small in this class of models is anomalous top decay, like t — >■ cr"'"r~. 
The other side of the coin is that if such decays are observed, the new physics must be something different from those 
described here, as the expectations will be in conflict with the i3-decay observables. 
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It is not yet certain whether there is a deficiency in the H — > r+r~ channel, but at the la level, the cross-section 
is slightly below the SM prediction. While it is too early to say anything about this channel, we would like to point 
out that the interactions discussed in this paper can potentially modify the predictions for this channel, without 
disturbing those for other channels. Further data from LHC will be eagerly anticipated. 
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